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Introduction to Density Functional Theory and 
Exchange-Correlation Energy Functionals 



R. O. Jones 

Institute for Solid State Research 
Forschungszentrum Jiilich 
52425 Jiilich, Germany 
E-mail: r. jones @fz-juelich. de 

Density functional calculations of cohesive and structural properties of molecules and solids 
can be performed with less computational effort than by using other methods of comparable 
accuracy. We outline the theory and discuss the local density approximations basic to most 
applications, and we discuss ways to modify these approximations to yield improved results. 

1 Foreword 

The NIC Winter School 2006 focuses on Computational Nanoscience, particularly on un- 
derstanding the properties of materials so that we can "tailor" new ones with novel prop- 
erties. A glance at the program shows the dominant role played by the density functional 
(DF) formalism in this area. I have been asked to present an introduction to this theory 
and, in a second lecture, to discuss approximations to the exchange and correlation energy. 
These lecture notes cover these topics, but the first lecture will follow a more historical (and 
more lively) approach to DF theory. The theory now plays a ubiquitous role in condensed 
matter calculations and in chemistry, but its general acceptance was slow in coming, and 
there has been much rewriting of history since its arrival. The following notes are derived 
in part from earlier articles 1 and provide a background for other talks in this School. 

2 Introduction 

The description of many-particle systems covers much of physics and chemistry. In com- 
putational nanoscience we focus on atoms, molecules, and solids, on their structural and 
cohesive properties, and on how they interact with each other. We usually emphasize the 
properties of specific materials rather than universal properties of all, and our interest on 
structural and binding means that the total energy E of a system and its components is 
very important. The density functional (DF) formalism is our main theme here, but we 
shall note how its perspective differs from that of other methods used to discuss many- 
particle systems. 

If we wish to discuss the properties of an interacting system, such as the elec- 
trons in a molecule or solid, it is natural to consider the many-electron wave function, 
^(^ljfbj —Tn)> where the r*j denote the particle coordinates and spins. This has been 
discussed in this School by Prof. V. Staemmler. Hartree 2 made one of the earliest approx- 
imations for \l> in writing the many-electron wave function as a product of single-particle 
functions, i.e. 

*(ri,r2,-) = i>i{ri) V'iv(rjv) (1) 
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Each of the functions tpi(fi) satisfies a one-electron Schrodinger equation with a potential 
term arising from the average field of the other electrons, i.e. 

(-^V 2 + Vext + $0 = eilfcO 8 ). ( 2 ) 

where the Coulomb potential $j is given by Poisson's equation, 

JV 

V 2 $ i = 47re 2 ^|V j | 2 , (3) 

3=1 

and T^xt is the potential due to the nuclei. Fermi statistics can be incorporated by replacing 
the product wave function by a single determinantal function. 3 4 This "Hartree-Fock" (HF) 
approximation leads to an additional, nonlocal exchange term in the Schrodinger equation, 
but does not change the single-particle picture, with the wave function described in terms of 
orbitals with particular spins and occupation numbers. The inclusion of "exchange" effects 
improves the total energy, and the HF approximation remains an indispensable benchmark 
in molecular physics. 

In 1960 Coulson, 5 however, noted that "it is now perfectly clear that a single configu- 
ration (Slater determinant) wave function must inevitably lead to a poor energy". A linear 
combination of such determinants would be better, and the inclusion of "correlation ef- 
fects" by this approach ("configuration interaction", CI) improves the many-particle wave 
function, from which most properties of interest can be calculated. Some such approaches 
will be discussed in this School by Dr. C. Hattig. The explosive increase in the number 
of configurations with increasing electron number means, however, that only systems with 
relatively few electrons can be calculated with high accuracy. The complexity of the re- 
sulting solutions can also hinder a simple interpretation of the results. 

Thomas 6 and Fermi 7 had proposed much earlier a scheme based on the density of 
electrons in the system, n(f). The Thomas-Fermi (TF) scheme assumes that the motions of 
the electrons are uncorrected, and that the corresponding kinetic energy can be described 
by a local approximation based on the results for free electrons, [n(r)]s . The potential 
of the system can be determined from solving Poisson's equation, and the requirement of 
constant chemical potential leads to the TF equation for n(f). This approach gives a poor 
description of real systems, but in its focus on the density it is a forerunner of later DF 
theories, and its mathematical properties have received much attention. 8 In general, the 
integrated quantity n(r) = N J dti- J df^^*(f, r 2 --f^)^(f, r 2 ...r*jv) should be easier 
to describe than the precise details of the wave function, \l>. 

Shortly after the original work, Dirac 9 proposed that exchange effects be included by 
incorporating a term derived from the exchange energy density in a homogeneous system, 

Dime 3 /3\* i 0.458 

£ Dirac _ _ _ n3 - a u (4) 

4 W r s 

where r s = (3/(47rn)) 1 /3 is a measure of the average interelectronic distance. The use of 
an approximate exchange potential in addition to the Hartree term in Eq. (2) was suggested 
by Slater, 10 who showed that the exchange potential in a system of variable density could 
be approximated by a term with the same density dependence. However, it is crucial 
to appreciate that the argument of Slater does not depend on the exchange energy in a 
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homogeneous system. The exchange interaction means that an electron of a given spin 
will be surrounded by a region where the density of electrons with the same spin (say, spin 
up) is reduced. The effect of the "exchange hole" can be estimated by replacing it by a 
uniform density inside a sphere of radius ro, and zero elsewhere. Since the hole contains 
a single electron, 10 ro = (3/(47rn-t-)) 1//3 , and we obtain an exchange energy with the same 
density dependence as that of Dirac. 

Gaspar 11 adopted the Dirac approximation to the exchange energy, and derived a 
slightly different exchange potential by using a variational approach. This method re- 
produced the HF eigenfunctions and eigenvalues for the Cu + ion very well, showing that 
an approximation based on results for a homogeneous system could give a satisfactory 
description even if the density is very far from uniform. 

The above work has been essential to the development of modern DF theory. As an 
approximation to the (nonlocal) exchange potential, the simple local density (LD) approx- 
imation provided the basis for much of the work on the electronic structure of solids in the 
1950's and early 1960's. The calculation of electronic band structures and related quan- 
tities is, of course, one of the success stories of modern condensed matter physics. The 
experience obtained in solving Schrodinger-like equations with local effective potentials 
and the improvement in numerical methods and computers during the past decades have 
been decisive in carrying out density functional calculations in recent years. 

Hohenberg and Kohn 12 showed that ground state properties, particularly the total en- 
ergy E, of a system of interacting particles could be related to the density distribution. For 
a nonmagnetic, nondegenerate ground state of a system of electrons in the presence of a 
scalar field V ext , they showed that the density uniquely determines the potential, up to an 
arbitrary additive constant. Their proof has been simplified and extended by Levy, 13 and 
the TF equation may be derived from this formalism as an approximation. Its application 
to the total energy of the interacting system results in a minimum principle, and the deter- 
mination of the energy of the ground state can be found by the solution of single-particle 
equations. 14 In practice, total energy calculations require approximations to be made for 
the exchange-correlation energy, E xc , and Kohn and Sham 14 showed that the local density 
(LD) approximation, 

£xc D = / dfn(^}e xc [n(^l (5) 

could be applied to the limiting cases of a slowly varying density and very high densities. 
Here £ xc [n] is the exchange and correlation energy per particle of a homogeneous electron 
gas with density n. The exchange contribution, e x , has the same form as that of Gaspar, 
and a variety of approximations exist for the correlation term, e c . 

The limiting cases noted above are not realized in atoms, molecules or solids, and Kohn 
and Sham commented that "we do not expect an accurate description of chemical bonding" 
with the LD approximation. 14 A decade passed before the first attempts were made to test 
its ability to describe the bonds in molecules, and it is remarkable that these showed that 
the LD approximation could generally reproduce ground state geometries, vibration and 
phonon frequencies, and moments of the density very well. The DF formalism remains 
the basis of most of the parameter-free calculations in extended systems, and has found 
widespread application to molecules and clusters. Most calculations have used the local 
density approximation [Eq. (5)] or modifications of it. 
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The practical necessity of approximating E xc leads to an essential difference in per- 
spective between the density functional and CI approaches. The goal of the latter is an 
exact numerical solution of the Schrodinger equation, which would yield exact answers 
for most quantities of interest, whereas even an exact solution of the density functional 
equations will only reflect the accuracy of the approximation for E xc . We shall discuss 
some of the many approximations that have been developed. In order to appreciate some 
qualitative features of the density in a real system, we show in Figure 1 the spherically 
averaged density in the ground state of the carbon atom. The density falls monotonically 
from the nucleus, with the radial oscillations often expected. The charge density in small 
molecules also tends to be relatively featureless, with maxima at the nuclei, saddle points 
along the bonds and a generally monotonic decay away from both. Furthermore, the elec- 
tron density in a molecule or solid shows relatively small departures from the overlapped 
densities of its atomic constituents. Energy differences result from subtle changes in rel- 
atively featureless density distributions, and it is not surprising that it is difficult to find 
simple prescriptions for determining them. 



C 2s V 





Figure 1. Spherical average of density in ground state of C atom 



3 Density Functional Formalism 

3.1 The Density as Basic Variable 

Ground state (GS) properties of a system, including the energy E, can be expressed as 
functionals of the GS electron density, i.e. they are determined by a knowledge of the 
density alone. We show that E[n] satisfies a variational principle. The Thomas-Fermi (TF) 
is a special case of this formalism. 

The basic theorems of the density functional formalism were derived by Hohenberg 
and Kohn, 12 and Levy 13 provided a simpler and more general derivation. We consider N 
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electrons moving in an external potential, V ex t(f), i.e. the Hamiltonian is 



N 



H = T + V ee + Y,Vext(T f i), (6) 



i=l 



where T and V ee are the kinetic and electron-electron interaction operators, respectively. 
For all "iV-representable" densities, n(r),° i.e. those obtainable from some antisymmetric 
wave function, ip(fi,r2, ■ ■■Fn), Levy defined the functional, 

F[n}= mm MT + VeeM, (7) 

ip— >n 

where the minimum is taken over all tp that give the density n. F[n] is universal in the 
sense that it refers neither to a specific system nor to the external potential V ex t (0 • If we 
denote E G s, tpcs, an d ^Gs(r) to be the ground state energy, wave function and density, 
respectively, then the two basic theorems of DF theory are: 

E[n] = J dfV^ t (f)n(f) + F[n] > E GS (8) 

for all 7V-representable n(r), and 

dfv(f)n GS (r) + F[n GS ] = E GS (9) 



To prove the variational principle [Eq. (8)] we introduce the notation V'min(r) f° r a 
wave function that minimizes F[n] in Eq. (7), so that 

F[n] = (^ in \T + V ee \^ in ) (10) 

Writing V = J2i V ex t(fi), we have 



/ 



drV^{¥)n{¥) + F[n] = (^ in \V + T + V ee \^ in ) > E GS , (1 1) 



according to the minimum property of the GS. This proves the inequality (8). Using the 
minimum property once more we find 

£ GS = {*1>GS\V +T + V ee |V>GS> < {^\V + T + VeeIVS (12) 

We subtract the interaction with the external potential and obtain 

<V>GS|T + V e e|V>Gs) < <«|T + V e e|V&. (13) 

On the other hand, the definition of V'mTn yields the reverse relation between the two sides 
of Eq. (13). This is possible only if 

<V> GS |T + Veelfe) = <CL S |T + Veeltt). (14) 



"Hohenberg and Kohn 12 worked in the space of V-representable densities, i.e. those that can be realized for 
some external potential, V e *t(r). This space is a subspace of iV-representable densities. 
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Then we have 

Egs = 

= J drV ext (f)n GS (f) + «9f|T + FeelO 

= J drV ext (f)n GS (7) + F[n GS ] (15) 

These are the basic theorems. An important result also follows from Eq. (14). If the 
ground state GS is nondegenerate, V'mfn = V'GS • If the ground state is degenerate, V'mfn 
is equal to one of the GS wave functions, and the others can also be obtained. The GS 
charge density then determines the GS wave function(s), from which all GS properties can 
be calculated. These properties are therefore functionals of the density. 

It has been said that Wilson 15 gave a shorter, if less general, proof of the definitive role 
of the electron density shortly after the DF formalism was first published. It is based on 
the "cusp condition" 16 satisfied by the density at the nucleus 

~d{n{r a ))~ 



Z - ~ l 

" ~ 2<n(0)> 



dr n 



(16) 



where the brackets (. . .) denote a spherical average. The density n then determines 
uniquely the charges Z a as well as the positions of the nuclei and hence the Hamiltonian. 

These theorems provide a general method for calculating GS properties. If an approx- 
imation can be found for F[n], we need to minimize E[n] in Eq. (13) for the potential 
v(r) of interest. This leads to corresponding approximations for Eqs an d ^Gs(r)- The 
same procedure applies to a ground state property X if we have an approximation for the 
functional X[n]. 

3.2 Thomas-Fermi (TF) Approximation 

Here the electrons are treated as independent particles, and the electron-electron interaction 
energy arises solely from the electrostatic energy, 

*e S [n] = £ / drj d? ^01. (17) 

Furthermore, we assume that the kinetic energy is given by 

T[n] = J drt[n(f)], (18) 

where t[n] is the kinetic energy density for a system of noninteracting electrons with den- 
sity n. We have 

i r ^ h 2 k 2 

(27r) 3 J\k\<k F 2w 
where 2(47r/3)fcjr 3 /(27r) 3 = n. This gives 

T [n] = C k J dr{n(f)}l, (20) 
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where C k = 3H 2 (3n 2 ) * /10m. 

We now minimize E(n) under the condition that the number of electrons is kept con- 
stant, 

j drn(r) = N. (21) 
Using the method of Lagrange multipliers, we require that 

E[n] +XN = T[n] + E es [n] + J dfn{f){V ext (?) + A} (22) 
has a minimum. The corresponding Euler equation is the TF equation, 

b -C k n{f) i + e 2 [ d? -^X- + V ext (f) + X = 0. (23) 
•j J \r — r'\ 

Decades of study have shown that the TF model provides a rough description of, for 
example, the charge density and the electrostatic potential, but it has severe deficiencies. 
The charge density is infinite at the nucleus, and it does not decay exponentially far from 
the nucleus of an atom, but as r~ 6 . It has also been shown that TF theory does not result 
in atoms binding to form molecules or solids. The absence of shell structure in the TF 
atom means that the observed periodic variation of many properties with changing atomic 
number cannot be reproduced. In fact, the atoms shrink with increasing atomic number Z 
(as Z~ 1 / 3 ), 11 and the model never predicts ferromagnetism. 18 



3.3 Single-particle Description of a Many-particle System 

Many of the drawbacks of the Thomas-Fermi approach can be traced to the approximate 
treatment of the kinetic energy. The task of finding good approximations to the energy 
functional is greatly simplified by using a different separation introduced by Kohn and 
Sham, 14 

E[n] = T [n] + J drn{¥) (v ext (7) + ^(f)) + E xc [n}. (24) 

T is the kinetic energy of a system with density n in the absence of electron-electron 
interactions, $ is the classical Coulomb potential for electrons, and the remainder E xc 
defines the exchange-correlation energy. T differs from the true kinetic energy T, but it 
is of comparable magnitude and is treated exactly in this approach. This removes many of 
the deficiencies of the Thomas-Fermi approximation, such as the lack of a shell structure 
of atoms or the absence of chemical bonding in molecules and solids. All terms in Eq. (24) 
but the exchange-correlation energy E xc can be evaluated exactly, so that the (unavoidable) 
approximations for E xc play a central role in the following discussion. 
The variational principle applied to Eq. (24) yields 
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where /j, is the Lagrange multiplier associated with the requirement of constant particle 
number. 6 If we compare this with the corresponding equation for a system with an effective 
potential V(r) but without electron-electron interactions, 



Sn(f) 6n(f) ~*~ ^' ^ 

we see that the mathematical problems are identical, provided that 

VW=7rt+4W + W (2?) 

Where the last term is referred to as the exchange-correlation potential V xc . The solution 
of Eqs. (26), (27) can be found by solving the Schrodinger equation for noninteracting 
particles, 

~V 2 + V(fij = 6*^(0. (28) 

yielding 

JV 

The condition (27) can be satisfied in a self-consistent procedure. 

The solution of this system of equations leads then to the energy and density of the 
lowest state, and all quantities derivable from them. The formalism has also been general- 
ized to the lowest state with a given set of quantum numbers. 19 In this case E xc depends 
on the values of the quantum numbers, and the density variations must remain within the 
space corresponding to the given quantum numbers. Instead of seeking these quantities by 
determining the wave function of the system of interacting electrons, the density functional 
method reduces the problem to the solution of a single-particle equation of Hartree form. 
In contrast to the HF potential, 

Vhf ip{r) = J dp Vhf(t, ?)*[>(?), (30) 

the effective potential, V(r) is local, and the equations are no more complicated to solve 
than Hartree's. The kinetic energy, the electrostatic interaction between core and valence 
electrons, and between valence electrons are treated exactly. Only the exchange energy, 
E x , and the even smaller correlation contribution require approximation. This is in marked 
contrast to the Thomas-Fermi and related methods, where the large kinetic energy term is 
approximated. 

We note here the problem of the "self-interaction correction" (SIC), where the Coulomb 
energy of interaction of an electron with itself must be cancelled by a contribution to E xc . 



b The functional derivative SF[n]/6n(f) is defined so that the first order change in the functional due to a change 
in its argument n(r) is: 



f SF 
5F = / dr——6n(r) 
J ontr) 
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The importance of SIC is obvious immediately if we consider a single-electron system 
such as the hydrogen atom. We also note the scaling condition on the exchange energy 20 

E x [X 3 n(\f]] = XE X [n(r)] , (31) 

which has proved to be useful in constructing functional approximations. 



3.4 Spin Systems 

The generalization of the single-particle equations to spin systems or an external magnetic 
field requires the consideration of the spin indices, a, of the one-electron operators ip a (f) . 
In general, this requires the replacement of V ext by a spin-dependent potential, V^(f), 
and the charge density n(f) by the density matrix 

n aP (f) = (9\^+(f)^ a (f)\9) (32) 

All ground state properties are functionals of the density matrix n a p, and E is stationary 
with respect to variations in n a p, provided 21 



dfn aa = N. (33) 



The application of the variational principle requires that n a p is 7V-representable, which is 
guaranteed for all non-negative, differentiable n a p that satisfy Eq. (33). 22 The effective 
spin-dependent potential in the single-particle equations is 

The original derivation of the Hohenberg-Kohn theorem leads to analytical problems 
in the case of spin systems. These authors showed that the ground state density n(f) deter- 
mined the ground state wave function \l/(fi , r*2, ...) uniquely, which determines the poten- 
tial V(f) up to an additive constant. For spin DF theory, however, von Barth and Hedin 21 
showed that the uniqueness of the potential is not guaranteed, and Capelle and Vignale 23 
found that this is true in general. A similar result was found by Eschrig and Pickett. 24 Kohn 
et al. 25 estimate that the practical consequences for spin DF calculations are not signifi- 
cant. Since the potential and magnetic field are unambiguous, and the kinetic, Coulomb, 
and exchange-correlation energies are defined in terms of but independently of n(r) 
and the spin density. A result that is useful in constructing functional approximations is 
"spin scaling", 26 which applies to the exchange energy in spin DF theory, 

£ X SDF H, n 1} = \ (£ x DF [2n t ] + £ x DF [2nJ) . (35) 

The approximation used most widely in total energy calculations of spin systems is the 
local spin density (LSD) approximation, 

£xc SD = / dfn(r) £xc [n t (f) , n x (r)} , (36) 

where £ xc [n^, n\\ is the exchange and correlation energy per particle of a homogeneous, 
spin-polarized electron gas with spin-up and spin-down densities and n^, respectively. 
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Numerous electron gas calculations have been performed over the years, and some parame- 
terizations of £ xc [n^, ri\\ will be discussed below. Different calculations lead to somewhat 
different results, but all are free of adjustable parameters. The Xa-approximation is ob- 
tained if £ xc in Eq. (5) is replaced by the exchange energy per electron and multiplied by 
3a/2, i.e. 

E^ a = -\aC J df {(n t (rO) 4/3 + (n^f 3 } , (37) 

1 /3 

where C = 3(3/47r) ' . The a-dependence of energy differences for a given atom or 
molecule is weak for values near 2/3, the exchange energy value (LSDX). There are sys- 
tematic differences between the two sets of calculations, with the latter approximation 
leading to overestimates of the relative stability of states with larger spin densities. 



3.5 Exchange-Correlation Energy, E xc 

The numerical advantages of the approach described are obvious. Efficient methods exist 
for solving single-particle Schrodinger equations with a local effective potential, and there 
is no restriction to small systems. We have noted, however, that E xc , is defined as the 
difference between the exact energy and other contributions that may be evaluated numeri- 
cally exactly. In practice, it is necessary to make approximations for this term, and we now 
examine it in detail. 

The crucial simplification in the density functional scheme is the relationship between 
the interacting system, whose energy and density we seek, and the fictitious, noninteracting 
system for which we solve Eqs. (28), (29). This can be studied by considering the interac- 
tion A/If — r'| and varying A from (noninteracting system) to 1 (physical system). This 
is done in the presence of an external potential, V\, 21 such that the ground state of the 
Hamiltonian 

H x = - 1 V 2 + Vext (0 +V X + AV ee (38) 

has density n(r) for all A. The exchange-correlation energy of the interacting system can 
then be expressed in terms of an integral over the coupling constant A. 19 ' 28 

E*c = \ [ drn(r) [ d? 1 - n^ c {r,r' - f), (39) 
* J J \r — r'\ 

with 

nxc (r, r'-f)=n(r')^ dX (g(f, r',X)-l). (40) 

The function g(r, r', A) is the pair-correlation function of the system with density n(f) and 
Coulomb interaction AV^, e - The exchange-correlation hole, n xc , describes the effect of the 
interelectronic repulsions, i.e. the fact that an electron present at the point f reduces the 
probability of finding one at r'. The xc-energy may then be viewed as the energy resulting 
from the interaction between an electron and its exchange-correlation hole. This picture is, 
of course, the same as that used by Slater. 10 
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We note three points: Firstly, since g(f, r') tends to unity as \ f — r'\ — > oo, the above 
separation into electrostatic and exchange-correlation energies can be viewed as an ap- 
proximate separation of the consequences of long- and short-range effects, respectively, of 
the Coulomb interaction. We may then expect that the total interaction energy will be less 
sensitive to changes in the density, since the long-range part can be calculated exactly. 

The second observation 19 arises from the isotropic nature of the Coulomb interaction, 
V ee , and has important consequences. A variable substitution R = r' — fin Eq. (39) yields 

E*c=lj dfn ^ I dR ^ ~R I dn nxc(r "' R) - (41) 

Eq. (41) shows that the xc-energy depends only on the spherical average of n xc (f,R), 
so that approximations for E xc can give an exact value, even if the description of the 
nonspherical parts of n xc is quite inaccurate. Thirdly, from the definition of the pair- 
correlation function, there is a sum-rule that requires that the xc-hole contains one electron, 
i.e. for all f, 



I 



dr' n xc (f,r' - r) = -1. (42) 



This means that we can consider — n xc (f, r' — r) as a normalized weight factor, and define 
locally the radius of the xc-hole, 

(43) 

\R\ 

This leads to 



<3>, = -/ 



= -\j dtn{f)( I) f . (44) 

showing that, provided the sum-rule [Eq. (42)] is satisfied, the exchange-correlation energy 
depends only weakly on the details of n xc . 19 In fact, we can say that it is determined by the 
first moment of a function whose second moment we know exactly. This argument simply 
makes more precise the ideas of Slater concerning the exchange hole and its consequences. 
The spin analogy of the exact expression for Ex C [Eqs. (39), (40)] is 

£ - = \ E / d?n a (r) [ dr'^^r f d\ (g a0 (f, r'; A) - 1). (45) 
a/3 J \r r \ Jo 

The atomic density distribution shown in Figure 1 makes it clear that arguments based 
on small departures from homogeneity cannot be applied. However, we have seen that a 
partial cancellation of errors can be expected provided that the sum-rule [Eq. (42)] is satis- 
fied. It is instructive to compare the form of the exchange hole found in LSD calculations 
with that determined exactly, i.e. in a Hartree-Fock calculation. This is shown in Figure 2 
for two representative values of f in the nitrogen atom. The approximate and exact holes 
are qualitatively different; the approximate hole is spherically symmetric and centred on 
f, while the exact hole is very asymmetric around its centre r'. The spherical averages 
are, however, remarkably similar, and the values of the exchange energy differ by only a 
few percent. The large differences in the exchange holes arise almost completely from the 
nonspherical components, which contribute nothing to the exchange energy. 
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Figure 2. Exact (solid line) and approximate (broken line) exchange hole n x (r, r' — r) [Eq. (40)] for a spin up 
electron in a nitrogen atom for r = 0.13 and 0.63 a.u. The top figure shows the hole along a line through the 
nucleus and the electron. The arrow indicates the nuclear position and f — r' = gives the electron position. The 
exact hole has a large weight at the nucleus, while the approximate hole is centered on the electron. The lower 
figure shows the spherical average of the hole around the electron. The area under the curve is proportional to the 
exchange energy. The figure also shows the value of defined in Eq. (43). 



4 Exchange- Correlation Energy Functionals 

The local density (LD) and local spin density (LSD) approximations have played a contin- 
uing role in DF calculations, and most functional approximations in current use are based 
on them. In this section we discuss these approximations as used in actual calculations, the 
reasons for the errors that result, and some of the modifications introduced to correct for 
them. 



4.1 LD and LSD Approximations 

The LD approximation [Eq. (5)] has an exchange component e x , given by the Dirac formula 
[Eq. (4)], and a much more complicated correlation component e c . The latter is commonly 
based on the results of quantum Monte Carlo (MC) calculations for a homogeneous elec- 
tron gas of different densities. 29 ' 30 The LSD approximation [Eq. (36)] also requires results 
for a spin-polarized system and a means of interpolating for one with partial polarization. 
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Perdew and Zunger 31 proposed the formula 



A In r s + B + Cr s In r s + Dr s r s <\, 



c "^/(l+Av^ + Ar,) r. >1. 

For dense electron systems (r s < 1) the random phase approximation provides the parame- 
ters for the LD approximation, and the parameters for a fully polarized gas can be obtained 
by scaling arguments. Other parameters were obtained by fitting to the results of Ceperley 
and Alder. 29 

Another form that has find widespread use is due to Vosko et al. 32 



£^ WN , ( r s \ 2b 

— = m —. ■ . H , arctan — — — - 

bx 



F(x ) 



In ( ^f^U 2(6 ~ 2X0) arctan ' 



F(y/K) J V4c - b 2 \V^ + b 



(47) 



Here F(x) = x 2 + bx + c and the fitting parameters, which differ for the polarized and 
unpolarized cases, are obtained from the data of Ceperley and Alder. 29 

Interpolation between unpolarized (U) and fully polarized (P) results is usually based 
on a formula involving the spin polarization (, 

Mn, C) = /(CK U C + (1 - /(C)) (48) 

where 

= n t (r) - H (f) 

Functional forms of f(() based on the Hartree-Fock and random phase approximations 
have been suggested by von Barth and Hedin 21 and Vosko et al., 32 respectively. 



4.2 Results of LD and LSD Calculations - General Trends 

We have noted above that the electron density in molecules and solids is generally close to 
a superposition of atomic densities. These are far from uniform and even have a cusp at the 
nucleus. The density is so far from uniform that the use of the LD and LSD approximations 
was treated initially by many with great reserve, particularly in the world of chemistry. We 
would not be discussing them today, however, if they did not have real value. In fact, they 
almost always lead to the correct picture of binding trends across and down the periodic 
table. Structures, bond lengths, and vibration frequencies in many systems are reproduced 
satisfactorily, and the deviations from measurements of these and other quantities are often 
quite systematic. We now discuss some of the deviations. 

Binding energies of molecules and cohesive energies of solids are usually overesti- 
mated. In weakly bound systems, e.g. where closed-shell molecules bind, these overesti- 
mates are often large, and the bond lengths are too short. In finite systems (atoms, surfaces, 
••• Wxc decays exponentially rather than with the inverse power of distance (the "image 
potential"). This is one of the consequences of the incomplete self-interaction correction. 
Negative ions (such as H~) are often unbound, although it may be possible to estimate 
binding energies in some cases. In semiconductors the energy gap between valence and 



13 



conduction bands is less than the measured values. This topic will be discussed in this 
school by Dr. A. Schindlmayr. The so-called strongly correlated systems, such as the in- 
sulating 3d-transition metal oxides, are not described well by LSD calculations, which 
predicts them to be metallic. Partially filled d- and /-bands, in general, require modifica- 
tions of the LD and LSD approximations. Such problems will be discussed in this school 
by Dr. S. Kurth. 33 Some calculated energy differences are in spectacular disagreement 
with experiment and provide us with a way to understand the defects of the LD and LSD 
approximations. We now focus on promotion energies in atoms, i.e. different states of 
lowest symmetry of the same atom. 



4.3 LSD Calculations - Sources of Error 

The above arguments show that the LSD approximation can give a reasonable description 
of systems where the density is far from homogeneous. As experience with molecular 
calculations developed, however, it became clear that LSD calculations gave rise to some 
persistent, occasionally systematic errors for a variety of systems. 1 The binding energies 
of sp-bonded molecules, for example, is often overestimated by ~1 eV per bond, and 
discrepancies of this size are simply unacceptable. 

When atoms bind to form a molecule or solid, there is generally a change in the relative 
s, p, d ... contributions to the wave function. Such energy transfers can be studied directly 
in atoms by calculating the energy required to change an s-electron into a p-electron, for 
example the 4 S (2s 2 2p 3 ) ground state of N into the 4 P (2s 1 2p 4 ) excited state. Since 
these are the lowest lying states of these symmetries, density functional calculations are 
justified. 19 

In Figure 3 we compare experimental sp-transfer energies for first row atoms and ions 
with the results of LSD calculations. 34 For atoms we define 

A sp = E(U 2 2s2p n - 1 ) - E(ls 2 2s 2 2p n - 2 ) (50) 

with a similar definition for the ions. While the overall trends are given well, the absolute 
deviations in AE sp are of variable quality, with particularly large deviations in the cases 
of O and N. Also shown are the HF 35,36 and LSDX (Xa) results. The trends are given 
fairly well by these calculations, but there are significant errors, and the comparison with 
experiment shows differences between the first and second halves of the row. To obtain 
insight into these problems we study some simple model systems, focusing on the exchange 
energy, where the presence of explicit formulas simplifies the discussion greatly. 
The HF exchange energy can be expressed in terms of exchange integrals, 37 



hj=et f if f d? ^m^mi t (51) 

J J I f — r' | 

where $j and $j are HF orbitals. Such an integral depends strongly on the nodal structure 
of and <£j. If and 3>j have different I- and m-quantum numbers, for example, 
the integrand oscillates and Iij is reduced. Since the LSD approximations for exchange 
(referred to as LSDX) and correlation are expressed only in terms of the charge density, 
we cannot expect them to include effects of the nodal structure in a precise way. 
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Be + B + C + N + + F + Ne + 



Figure 3. The sp transfer energies A sp for the first row atoms and ions, (a) experimental and LSD results, (b) 
Hartree-Fock and LSDX (Xce) results. The energies are in eV. 



As an example, we consider the sp-transfer in an F atom, where an s± electron is 
transferred to a p± orbital, 34 

u 2 2 S (n)2 P (mu) -> u 2 2 S tf)2p(mm) ■ (52) 

The change in the exchange energy due to this process is 37 



(53) 



where the Slater integrals G k are defined by 

P OO P OO }x 

G k (i,j)=e 2 J o drr 2 J o dr'(r') 2 ^ 



4>i{r)4>j{r)4>i{r')4>Ar') ■ (54) 



Here r< (r>) is the smaller (larger) of r and r', and 4>i(r) is the radial part of 3>j(f). 
The first term in Eq. (53) is the exchange interaction between a p-electron with m = — 1 
and two p-electrons with m = 1 and 0. The second term is the interaction between an s- 
electron and two p-electrons. The integrand of Eq. (51) corresponding to the p-p interaction 
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has two nodal planes as a function of r, while there is only a single nodal plane for the s-p 
interaction. The latter is larger as a result. Using realistic values for the Slater integrals, 
we find that AE X ~ 6 eV. The LSDX calculation, however, gives similar radial extents for 
the s- and p-orbitals. If we assume that the radial parts are identical and neglect the small 
nonspherical corrections, the LSDX approximation predicts that the exchange energy is 
unchanged by the sp-transfer, and the LSDX value for A sp differs from the HF result by 
6 eV. The LSDX and LSD (both exchange and correlation) approximations show similar 
trends, although the deviation between the LSD result and experiment is reduced to 2.6 eV, 
since the large change in the exchange energy is compensated in part by a change in AE C 
of opposite sign. 

A simple model problem. In order to gain insight into these problems, Gunnarsson 
and Jones 34 performed extensive LSDX and LSD calculations for atoms and compared 
the results with HF calculations and experiment, respectively. The differences between 
the LSDX and HF results are reproduced remarkably well by model calculations where 
we isolate the effects of the /-dependence of the orbitals by assuming that the s, p and d 
orbitals have the same radial dependence. With this assumption G k only depends on k. 
For the 3s orbital in calcium, the numerical relations 

G o = _G^ = _G^ = _G^ = _&_ (55) 
0.680 0.516 0.414 0.344 

and 

' df e x [n{f), 0] =0.451G° (56) 



were found, where n(f) is the charge density due to a Ca 3s electron. Similar coefficients 
were obtained for the orbitals of other elements of the iron series. Eq. (56) shows that the 
unphysical self-interaction, —G°/2, is cancelled to within about 10%. In Figure 4(a) we 
compare the LSDX and HF exchange energies as a function of the number N of spin up 
electrons. In the HF case, the shells are filled in the order s, p, d and within each shell 
in the order m,m — 1, .., —m. In the LSDX case the small nonspherical contributions are 
neglected. The LSDX and HF results agree remarkably well. 

Figure 4(b) shows results for the interelectronic exchange energy, for which the self- 
interaction has been subtracted. With the filling order s, p, d, the LSDX and HF results are 
in rather good agreement, with the magnitude of the LSDX results being somewhat larger. 
Figure 4(b) also shows results for occupations where a subshell is left empty. This does 
not influence the LSDX results in the present model, but it has a pronounced effect on the 
HF results. For instance, the curve "s-shell empty" is lower than the curve corresponding 
to the s, p, d filling, since the p-p exchange interaction is smaller than the s-p interaction. 
It follows from the orthogonality of the HF orbitals that the exchange hole contains one 
electron, and a similar sum rule is also satisfied by the LSD and LSDX approximations. 19 
Aspects of orthogonality and node formation are then included in all these schemes, al- 
though Figure 4(b) shows that the sum rule does not guarantee a good description. 

This model calculation illustrates two important conclusions: (i) If we occupy the or- 
bitals with the minimum number of nodal planes consistent with the Pauli principle, the 
trends in the interelectronic exchange energies are reproduced well by the LSDX approxi- 
mation. The absolute value is overestimated in all systems considered by Gunnarsson and 
Jones. 34 (ii) The energy for the transfer from such a state to a state with one additional 
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Figure 4. Exchange energies in model described in text, (a) SD (dashed curve) and E™ (solid curve) as a 
function of the number of electrons N (in units of — G°); (b) the interelectronic exchange iS x> ; nt for these two 
approximations. In the HF approximation, we show the dependence for different schemes of occupying the 
orbitals. 



node is often underestimated substantially in the LSDX approximation. These results were 
derived from atomic calculations, but the arguments can be extended to small molecules. 34 
This shows that the nodal structure of the wave functions can greatly affect the accuracy 
of the LSD approximation, and we have identified classes of problems where the LSD 
results must be treated with caution. For states with the minimum number of nodal planes 
consistent with the orthogonality of the orbitals, the LSD approximation usually leads to 
a moderate overestimate of the exchange-correlation energy. For states with additional 
nodal planes the exchange-correlation energy is often greatly overestimated. In atoms, the 
depopulation of s-orbitals can lead to large errors, and similar effects may be expected in 
bonding situations where sp or sd hybridization reduces the s occupancy. 



4.4 A Simplification of LD Calculations: Harris Functional 

Before discussing some of the suggestions for overcoming problem with the LD and LSD 
approximations, we shall mention one way of simplifying LD calculations. In 1972, Kim 
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and Gordon 38 39 suggested the following extension of the Thomas-Fermi approach: The 
density of a system comprising closed-shell atoms or molecules is approximated by over- 
lapping the HF densities of the free atoms (or molecules), and the energy is then calculated 
using the TF approximations for the electrostatic [Eq. (17)] and kinetic energy [Eq. (18)] 
terms and Eq. (5) for the exchange-correlation energy. The energies of the constituents 
are estimated in the same way, and the binding energies of rare gas molecules so obtained 
agree remarkably well with experiment. 

The same approach has been studied in the DF context by Harris. 40 For systems where 
the density is close to the overlapped density of the constituents, one can derive an expres- 
sion for the binding energy for a given geometry, 



Er ~ a n~*n ~ J df nj (0 l -4>f (0 + V% (r) 



+ E xc [n f ]+E§, (57) 



where rif is the sum of the densities of the fragments, <j>f and V£? are the corresponding 
Coulomb and exchange-correlation potentials, and En is the internuclear repulsion. The 
e n are eigenvalues of the potential 

V{¥) = </> f (f) + V x 7 (f) + Vext (0 , (58) 

where V ex t is the external potential. This expression has the same form as that of Kim 
and Gordon, with the important difference that the change in kinetic energy (for bonding 
between atoms A and B with potentials Va and Vg) is 

6T, = 2^<5g„ - f dfn(n) [V{¥) - V A (n)} 

- J drn{r 2 )[V{¥) - V B {f 2 )], (59) 

instead of the corresponding Thomas-Fermi expression. 

The results of this scheme agree satisfactorily with self-consistent DF calculations for 
He2, Be2, C2, N2, F2 and CU2. 40 It provides a useful alternative in cases where more 
detailed calculations would be prohibitive, and it has been adopted in some molecular 
dynamics calculations, where more detailed calculations are performed at regular intervals 
to ensure that the trajectories are reliable. It is, however, not self-consistent, and spin plays 
no role. 



4.5 Modifications of the LSD Approximation: Pair-correlation Functions 

We have seen that the total energy of a system is insensitive to many details of the pair- 
correlation function g(f,r'), which occurs in Eq. (40). It is nevertheless tempting to 
modify the pair-correlation function to remove some of its more unpleasant features. The 
expression for the exchange-correlation energy [Eq. (41)] shows that the differences be- 
tween the contributions of the exact and approximate exchange holes are largely due to 
the non-spherical components of the hole. Since these do not contribute to E xc , total en- 
ergies and total energy differences can be remarkably good, even in systems where the 
density distribution is far from uniform. In the LD approximation, we assume that the xc- 
hole n xc (r, r — r') depends only on the charge density at the electron. It would be more 
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appropriate to assume 41 that n xc depends on a suitable average n(f), 

n xc (r,r' - r) = n{r) I d\ [gh(r- r',\,n(r)) - 1]. (60) 



It is possible to choose the weight function which determines n(r) so that the functional 
reduces to the exact result in the limit of almost constant density. The approximation 
[Eq. (60)] satisfies the sum rule [Eq. (42)]. Other prescriptions for the weight function 
have been proposed. 

An alternative approximation is obtained if we keep the proper prefactor n(r') in 
Eq. (40), leading to the so-called "weighted density" (WD) approximation, 

n xc (f, r' — r) = n(r') G(\f — r')\,n(r)), (61) 

where n(f) is chosen to satisfy the sum rule [Eq. (42)]. The WD models provide a link 
between e xc , the density-density correlation function, and the response function of a homo- 
geneous system. As an example, we now discuss an analytic form 42 that is computationally 
simpler and can be chosen to give exact results in certain limits. We assume that 

G(f, n) = C(n) (1 - exp[-A(n)/|r| 5 ]) , (62) 

where C and A are parameters to be determined. The functional G behaves as |r|~ 5 for 
large distances, which is needed to obtain an image potential. For a homogeneous system 
with density n, we require that the model functional should both fulfil the sum rule for 
n(f) = n and give the exact xc-energy. This leads to two equations 

J dfG(\r\,n) = -l (63) 

y / df^G(\r\,n)=e xc (n), (64) 

which are sufficient to determine the two parameters C(n) and A(n). 

This functional is exact in several limiting cases: (i) for a homogeneous system, (ii) 
for one-electron systems such as the hydrogen atom it gives an exact cancellation of the 
electron self-interaction, (iii) for an atom it gives the correct behaviour of the xc-energy 
density far from the nucleus. £ xc (r) = — e 2 /2r, (iv) far outside the surface it gives the 
image potential e xc {z) = — e 2 / (4z). The LSD approximation gives qualitatively incorrect 
answers for (iii) and (iv), and the cancellation in (ii) is only approximate. Since (ii) is 
satisfied, this approximation provides a "self-interaction correction" in the sense we discuss 
below. The errors in the WD approximation for the total energy of atoms are about an order 
of magnitude smaller than those in the LSD approximation, but the sp- and sd-transfer 
energies are little improved over LSD results. 1 

4.6 Self-interaction Corrected (SIC) Approximations 

In the DF formalism each electron interacts with itself via the Coulomb electrostatic en- 
ergy, and this unphysical interaction would be exactly cancelled by a contribution from 
the xc-energy. In the LSD approximation this cancellation is imperfect, but numerically 
rather good. The incorrect treatment of the self-interaction in approximate functionals has 
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led some workers to consider self-interaction corrected (SIC) functionals. Within the LSD 
approximation the SIC functional takes the form 

Esic = £lsdK (r) ; "4.(r)] - X!^' (65) 

ia 

where Elsd is the LSD energy functional and 

5 ia = e l J drj dr' " ig( ^: (r,) + J^fW 0) (66) 

is the self-interaction correction (SIC) for the orbital i with spin a and charge density 
n-ia (r) ■ The first term in Eq. (66) is the self-interaction energy, and the second is the LSD 
approximation to the xc-energy of a fully spin-polarized system with density rii a (f). This 
functional is exact for a one-electron system. It leads to the equation 

[- ^V 2 + V(r) + V* c (r)]<f> ia = J] Kj4>ja , (67) 

i 

where V(f) is the effective potential (27) entering in a normal LSD calculation, and 
V^ c (f) is additional potential resulting from the term Eq. (66). This potential is orbital- 
dependent and the Lagrange parameters A^- are introduced to ensure that the solutions 
4>i (r) (referred to as "local orbitals") are orthogonal. The SIC approximation is not invari- 
ant under a unitary transformation of the orbitals, and different basis sets lead to different 
total energies. 

Unfortunately, ionization and transfer energies of atoms are not generally better in the 
LSD-SIC approximation, at least if non-spherical corrections are neglected. The improved 
total energies in SIC calculations can, in fact, be traced to a much better treatment of 
the innermost core electrons, 43 which play a relatively minor role in most chemical and 
physical processes of interest. 



4.7 Modifications of the LSD Approximation: Gradient Expansions 

Hohenberg and Kohn 12 introduced the local density approximation in the context of DF 
theory, but they also pointed out the need for modifications in systems where the density is 
not homogeneous. One such modification was the approximation 




where the kernel K xc is related to the dielectric function of a homogeneous medium. This 
approximation is exact in the limit of weak density variations 

n(r) =n + An(r), (69) 

where | An(r) | <C no, but the results for real systems were not encouraging. For free atoms 
the energy is infinite, 41 indicating that the sum-rule [Eq. (42)] is not satisfied. In spite of 
this, gradient expansions have played an extremely important role in DF theory and its 
applications. We shall now discuss some of the forms used. 
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4.7.1 Generalized Gradient Approximations (GGA) 



The second-order generalized gradient approximation to the exchange-correlation energy 
is written 

Eg GA (n) = J dre xc (n, |Vn|, V 2 n) . (70) 

Many forms of such functionals have been suggested, and some of the most important (and 
all their parameters) are given by Filippi et al. 44 These authors also compared these forms 
with the exact results for an exactly soluble system of two electrons in an external harmonic 
potential. We shall discuss some functionals that are widely used in actual calculations. We 
use some standard abbreviations: 



For the exchange contribution, Becke 45 suggested the following form: 

P x 2 



£ -£ LD 

tx — t„ 



1 



(72) 



2 1 / 3 A X 1 + Q/3x arcsinh(x) J ' 

where x = 2(6tt 2 ) 1/3 s = 2 1/3 |Vn|/n 4/3 , A x = §(3/7r) 1/3 , and j3 = 0.0042. The 
parameter j3 was optimized to give exchange energies of noble gas atoms using Hartree- 
Fock orbitals). This functional has been coupled with the LYP functional for correlation: 46 

{ n + bn ^ 2/3 Cfti 5/3 - 2t w + ^ (t w + \v 2 n) exp(-cn- 1/3 )J 

(73) 



>n =lf--V 2 « (74) 



where 



n 

and C F = 3/10(3tt 2 ) 2 / 3 , a = 0.04918, b = 0.132, c = 0.2533, d = 0.349. It surprised 
many that the combination of Becke exchange with either the VWN or LYP correlation 
functionals outperformed correlated ab initio methods (MP2 and QCSID) in calculations 
of atomization energies for 32 molecules. 47 The importance of using a non-local exchange 
has been underscored by many subsequent calculations 48 

The question of ways to improve exchange-correlation functionals has led to a clear 
dichotomy of views. The first implies that the search for a functional is so difficult that we 
should develop a reasonable form and fit its adjustable parameters to experimental data. 
The BLYP and B3LYP functionals are of this type. From this standpoint, DF theory is then 
"of semi-empirical nature". 49 Experimental data from up to 407 atomic and molecular 
systems have been used, 49 in this case with a functional with 15 adjustable parameters. 
The second approach seeks to build on the positive features of the LSD approximation 
by incorporating exact constraints and hoping that the added features will automatically 
improve the description of reality. We shall focus on functionals in the second category. 

Perdew, Burke, and Ernzerhof 50 have described a functional form (PBE) that has sev- 
eral attractive features: 

E xc = J drn{f)e^ (n(f))F xc (n, (, s), (75) 
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where 



F x (s) =1 + K- 



H = /3(n 2 /3) = 0.21951, and /3 = 0.066725. The correlation energy has the form 



E c = J drn{f) [^ SD (n,C) +H(n, (,t)] 



with 



H(n, (,t) = 



1 + 



2 r 



1 + At 2 



l+At 2 + AH 4 



]} 



(76) 



(77) 



(78) 



The spin-scaling factor 0(C) = [(1 + C) 2/3 + (1 - C) 2/3 ] A 7 = (1 
0.031091, and 

A = £ [exp (- e L D [n]/( 7 3 e 2 /«o)) - if' ■ 



ln2)/7r 



2 _ 



(79) 



The PBE form has several advantages, including: (1) In the slowly varying limit 
(t — > 0), H is given by the second-order gradient expansion, (2) In the rapidly varying 
limit (t — > oo), H — > —£^ D , and correlation vanishes, (3) It cancels the logarithmic singu- 
larity of e[r D in the high density limit. The accuracy of the PBE functional for atoms and 
molecules has been compared with results of LSD, BLYP, and B3LYP functions has been 
made by Ernzerhof and Scuseria. 51 The PBE functional performed as well as B3LYP for 
the properties considered by these authors. 
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Figure 5. Left: Structure of crystalline polyethylene: (a) normal to chain axis, (b) along chain axis. Right: 
Structural parameters of crystalline PE from x-ray diffraction, and calculated structures for a single chain. Bond 
lengths are in A, angles in degrees. Values marked with an asterisk (* ) were assumed, and not measured directly. 
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Generalized gradient approximations generally lead to improved bond angles, lengths, 
and energies. In particular, the strengths of hydrogen bonds and other weak bonds between 
closed shell systems are significantly better than local density results. However, the self- 
interaction problem remains, and some asymptotic requirements for isolated atoms are not 
satisfied (V® GA vanishes exponentially far from the nucleus, V GGA (f — > 0) — ► — oo). 

An example of the differences between LSD and PBE calculations is shown in 
Figure 5. 55 Polyethylene crystallizes in a structure with parallel chains, and the table shows 
that the structure of the individual chains is described rather well, with the C-C bonds being 
slightly longer in PBE and BP 56 calculations than in those using LD. There are, however, 
striking differences in the interchain separations. The LD calculations lead to a pronounced 
minimum with a shortest distance between C atoms in different chains (3.67 A) less than 
the measured value (4.59 A), 54 while PBE leads to a modest overestimate (4.66 A). 

4.7.2 Meta-GGA 

The next step in the development of gradient approximations is to incorporate the kinetic 
energy density (or the Laplacian of the density). A version based on the PBE form was 
described by Perdew et al., 57 who added the variables r-j- and rj., where 

occ. 1 

^(^L^w^i 2 (80) 

i 

is the kinetic energy density for the occupied Kohn-Sham orbitals. However, this and other 
forms initially developed included parameter(s) found by fitting to experimental data. This 
last feature was avoided in the recent work of TPSS, 58 whose form satisfied the requirement 
that the exchange potential be finite at the nucleus for ground state one- and two-electron 
densities. This is a constraint satisfied by LSD, but lost in GGA. Extensive numerical tests 
for atoms, molecules, solids, and jellium surfaces showed generally very good results. 

4.7.3 Hybrid Schemes: Combination of HF and DF Schemes 

In section (4.3) we saw just how poor the exchange energy differences could be between 
states whose wave functions have different nodal structures. It has also been noted for 
many years (see, for example, Ref. 27) that errors in the LD descriptions of exchange and 
correlation tend to balance. This suggests that a combination of Hartree-Fock exchange 
and DF calculations could be useful. 

E% hM = aE™ + E c , (81) 

where a can be chosen to satisfy particular criteria. A formal justification for such hybrid 
schemes was given by Gorling and Levy. 59 The B3LYP functional, 48 which is widely used 
in the chemical community, is an example. Here three adjustable parameters are used to fit 
calculated values to a molecular data base. 

4.7.4 Optimized Effective Potential (OEP) 

In the Kohn-Sham implementation of DF theory we describe the kinetic energy functional 
for the system of non-interacting electrons in terms of the orbitals ip ia (f) . The best local 
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effective potential (a requirement of DF theory) can be determined variationally for each 
a such that, when substituted into a stationary single-particle equation it leads to a set of 
eigenfunctions that minimize E[ipi a (r)]. This was first suggested by Sharp and Horton 60 
following the local density description of exchange, 10 and it is referred to as the optimized 
effective potential (OEP). 61 Its determination requires the solution of an integral equation 
for each V XC:a (r), which is a demanding exercise. How this can be done, as well as possible 
simplifications, are discussed in this School by Dr. S. Kurth. 

5 Concluding Remarks 

The density functional theory is the basis of most of the calculations in computational 
nanoscience that are free of adjustable parameters. It is clear from the present discus- 
sion that "DF calculations" cover a colourful variety of functional approximations, and 
not all are "free of adjustable parameters". The development of approximations to the 
exchange-correlation functionals over the past 20 years has improved the performance of 
DF calculations, and John Perdew is no doubt not alone in thinking that progress up the 
"Jacob's ladder" 62 (Figure 6) will continue until energy differences can be determined to 
within ~ 1 kcal/mol ("chemical accuracy"). The numerical cost increases as one climbs, 
and this may not necessarily bring more information. Most problems in "computational 
nanoscience" are performed part way up the ladder, and this situation will probably remain 
true for some time. 

DF calculations should always be performed with all critical facilities intact. The fa- 
miliar program "packages" will almost always lead to an answer, but it is easy to obtain 
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Figure 6. "Jacob's ladder" of DF schemes according to Perdew and collaborators. 
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answers void of physical content. It is odd to have to conclude this manuscript on such a 
note, but years of refereeing articles in the field suggest that it is necessary. 
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